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The microscopic structure of spin-orbit fields for the technologically important Fe/GaAs interface 
is uncovered from first principles. A symmetry based method allows to obtain the spin-orbit fields — 
both their magnitude and orientation — for a generic Bloch state, from the electronic band structure 
for any in-plane magnetization orientation. It is demonstrated that the spin-orbit fields depend not 
only on the electric field across the interface, but also surprisingly strongly on the Fe magnetization 
orientation, opening prospects for their magnetic control. These results give important clues in 
searching for spin-orbit transport and optical phenomena in ferromagnetic/nonmagnetic systems. 
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In solid-state systems lacking space inversion sym- 
metry spin-orbit coupling (SOC) acts on the electronic 
structure as a spin-orbit field (SOF), which is an effec- 
tive magnetic field whose direction and magnitude de- 
pend on the electron momentum fP, ^ . The most promi- 
nent examples are the Dresselhaus spin-orbit field [3 de- 
scribing the effects of bulk inversion asymmetry (BIA) 
in zinc-blende semiconductors, and the Bychkov-Rashba 
spin-orbit field [4 , describing the effects of structure in- 
version asymmetry (SI A) in asymmetric quantum wells. 
Apart from semiconductor structures, where Bychkov- 
Rashba coupling has been extensively studied 21 13 E] 
has been investigated in many other systems, for example 
on metallic surfaces [7]4T2^. graphene on a Ni substrate 
[13], or in Au and Ag monolayers on W(llO) substrates 
[14]. A striking manifestation of spin-orbit coupling in 
condensed matter is the spin-momentum locking in topo- 
logical insulators . 

Spin-orbit coupling can be controlled by an electric 
field [T6 . This fact has for long been used to moti- 
vate spintronics applications as epitomized by the Datta- 
Das transistor [17] in which the gate controls the spin- 
orbit induced spin precession of the itinerant electrons 
in a transistor channel. But spin-orbit coupling is also 
important for anisotropic magnetotransport. Tunneling 
anisotropic magnetoresistance (TAMR), for example, can 
be used to control electrical transport by rotating the 
magnetization orientation of a single ferromagnetic layer. 
It has been observed and studied in a variety of sys- 
tems, GaMnAs/Al [18], Fe/GaAs, [3 US], CoFe/GaAs 
[20] (inserting an MgO barrier suppresses TAMR here 
[21], a clear evidence for interface induced symmetry of 
the effect), Co/Ft [22], Si/ferromagnet junctions [23], 
resonant tunnel devices [24], or on an atomic scale in 
STM experiments [25]. Interfacial spin-orbit coupling 
has been proposed to control thermoelectric anisotropics 
in helimagnetic tunnel junctions [26] and produce spin- 
transfer torque in ferromagnet-topological insulator junc- 
tions [27]. 

In earlier studies of spin-orbit coupling on surfaces [7]- 
HHIKIIEH] and interfaces [Tl[14] the spin-orbit Hamil- 



tonian was extracted by fitting the energy bands close 
to the F point assuming a Bychkov-Rashba-type cou- 
pling. This standard procedure requires a priori knowl- 
edge of the specific functional form of the spin-orbit field 
and applies only to very small k- vectors for which small- 
momentum expansions are meaningful. Here we intro- 
duce a novel method to obtain spin-orbit fields (not just 
the functional parameters) for a generic k-point directly 
from ab-inito data. On the example of an Fe/GaAs junc- 
tion, important for room temperature spin injection [29]- 
i33j and TAMR [191 El] : we derive a formula for the mag- 
nitude and direction of the momentum dependent spin- 
orbit fields directly form the electronic band structure. 
The results show highly anisotropic (with respect to the 
momentum orientation) patterns, which take on different 
forms, from the ones known in semiconductor physics for 
small momenta to more exotic ones for Bloch states fur- 
ther away from the F point. 

One fascinating outcome is a qualitative dependence 
of the spin-orbit fields patterns on the band (energy), 
consistent with the bias-induced inversion of the TAMR 
observed in experiments [19l[35|. Even more important, 
in addition to their sensitivity on an electric field, the 
spin-orbit fields can depend unusually strongly on the 
magnetization direction, to the point that the anisotropy 
axes can be flipped by rotating the magnetization. We 
emphasize that those effects are caused by the symmetry 
of the interface, not of the bulk structures, making them 
particularly important for lateral transport in ultrathin 
hybrid ferromagnet-nonmagnet junctions. 

We consider thin Fe/GaAs slabs. The small lattice 
mismatch between Fe (2.87 A) and GaAs (5.65 A) allows 
for a smooth epitaxial growth of Fe on a GaAs (001) 
surface. Early investigations of the stability of 1 x 1 
Fe/GaAs interfaces within density functional theory [36] 
showed that when more than two atomic layers of Fe are 
deposited on a GaAs (001) surface, the flat or partially 
intermixed interfaces are more stable than the fully inter- 
mixed one, the As-terminated flat interface being more 
stable that the partially intermixed one. On the other 
hand, a recent Z-contrast scanning transmission electron 
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microscopy reported a single plane of alternating Fe and 
As atoms at an Fe/AlGaAs interface [371 [38]. Since the 
choice of the interface is not important to the message of 
our paper, we choose an As- terminated fiat interface. 

The electronic structure of an ideal Fe/GaAs slab, con- 
taining 9 (001) atomic layers of GaAs with the diagonal 
lattice spacing d = a/ a/2 = 3.997 A and three atomic 
planes of bcc Fe, has been calculated using the full po- 
tential linearized augmented plane wave technique imple- 
mented in the FLEUR code [39 and a generalized gra- 
dient approximation for the exchange-correlation func- 
tional [40]. The SOC has been treated within the second 
variational method. 

The band structure of the Fe/GaAs slab along the high 
symmetry lines connecting the S — F — X points in the 
Brillouin zone (BZ) is shown in Fig. [l] for a magnetiza- 
tion orientation along the [110] direction. The spin char- 
acter of bands 1 and 2 in Fig. [l] is basically determined 
by the interface atoms. The interface unit cell contains 
interfacial As, the neighboring Ga, and two Fe atoms. 
The spin-up character of band n = 2 is dominated by 
the interfacial As atom, its neighboring Ga atom and Fe 
atom above Ga, while the spin-down character of band 
n = 1 comes mostly from the two Fe atoms. 



around the F point as 




FIG. 1: Calculated band structure for the Fe/GaAs slab and 
magnetization along [110]. The states with spin-up (spin- 
down) character at the Fe/GaAs interface are emphasized by 
blue filled (red open) circles whose radii are proportional to 
the corresponding charge density at the interface atoms. The 
inset shows the As-terminated flat 1x1 interface model as- 
sumed in the study. 

The non-centrosymmetric GaAs layer is of D2d sym- 
metry, exhibiting the BIA spin-orbit coupling. The in- 
terface lowers the symmetry to C2v with the twofold ro- 
tation axis C2 along the growth direction [001] [2 . The 
C2V symmetry accounts for both the BIA and SIA; the 
C2V spin-orbit field lies in the plane of the slab, perpen- 
dicular to the growth direction. Since C2V symmetry has 
only one-dimensional orbital irreducible representations, 
away from accidental level (anti) crossings the spin-orbit 
fields (even at high symmetry points) can be described 
by spin 1/2 Pauli matrices. 

The most general SOC Hamiltonian consistent with 
C2V symmetry can be written for the in-plane momenta 



0)kycrx , (1) 



where kx and ky are the components of the in-plane wave 
vector k, cFx and cFy are the Pauli matrices, and x and 
y correspond to the diagonal [110] and [110] crystallo- 
graphic directions in GaAs, respectively; refers to the 
magnetization direction with respect to the [110] crys- 
tallographic direction of GaAs and n labels the band of 
interest. The functional parameters fin and rjn^ 
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are even in the momenta and, what is crucial and new 
here, depend in general on the magnetization direction. 

The values of the expansion parameters /in\ r]n^ {i = 
0, 1, 2, . . .) determine the specific form of the SOF. For ex- 
ample, if /J.l^^ an and r^i^^ -a^ (/in ^ r]n-^ Pn), 
Hso reduces in the limit of small k = |k| to the well 
known Bychkov-Rashba [4 (linearized Dresselhaus [3]) 
SOC with an (Pn) denoting the Bychkov-Rashba (Dres- 
selhaus) SOC parameter of the nth band. By introducing 
the SOF field 
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Eq. ([T]) can be rewritten as Tiso = '^n(k) • tr, where a is 
the vector of the Pauli matrices. 

We first analyze the spin-orbit fields for in-plane mag- 
netization directions. Since the exchange field dominates 
over spin-orbit coupling and the magnetization lies in the 
plane of the layers, the SOC contribution to the energy 
can be treated within first order perturbation theory. 
From the symmetry properties we find (see the Supple- 
mentary Material for the details) the following relations. 



,(k,^) 



A^-(k,^) + r-(k,( 



2 cos (9 



and 



where 



(k,^) 



AEf,°{k,0)-Tl-{k,0) 
2 sin 61 



E„ik,0)-En(--k,0) 



(4) 



(5) 



(6) 



r-(k,( 



kxjkyjO) En(kx-) ky-)^ 



(7) 



and (T refers to the spin character of the nth band. The 
above relations allow us to extract the components of 



the SOF directly from the ab-initio energy bands. In the 
particular cases of ^ ^ 7r/2 and 6^0 the numerators and 
denominators in Eqs. Q and ([5|, respectively, vanish. In 
such cases the SOF is obtained by L'Hopital's rule. The 
validity of Eqs. Q and ([5| is not restricted to the vicinity 
of the F point but holds also for large momenta. The only 
restriction is that the /c-space region of interest must be 
away from energy anticrossings. 

Figure |2] establishes the proof of principle for the mag- 
netization dependence of SOFs. It shows the SOF, 
ii;(k) (bottom parts), and polar plots of its strength 
w = |i^(k)| (upper parts), for the interface band n = 1. 
The SOF is computed on three different contours around 
the F point, k = 7r/100(i, ir/Sd and n/bd and plotted 
in Figs. [2|i), b), and c), respectively. The left (right) 
panel corresponds to the magnetization pointing along 
[110] ([110]). The C2v symmetry of the SOF is pre- 
served for all k. In particular, close to the F point 
the SOFs resemble the interference of Bychkov-Rashba- 
type and Dresselhaus-type SOCs (see Fig. [2^). However, 
away from the F point higher in k terms become relevant 
and more exotic patterns — we call them spin-orbit-field 
butterflies — in the SOF appear (see Fig. [2]3, c). The lin- 
ear terms are dominant up to about 5% from the BZ cen- 
ter, where the SOF exhibits a very strong dependence on 
the magnetization orientation. Note that the principal 
symmetry axes of the SOF can even be flipped by turn- 
ing the magnetization orientation. This remarkable effect 
opens the perspective of a magnetic control of spin-orbit 
fields. 

Close to the F point the SOF is determined by the 
contributions linear in the wave vector components kx 
and ky and characterized by Bychkov-Rashba-type and 
Dresselhaus-type SOC parameters, = [fin^ — ^^]/2 
and /3n = [fin'' +r/i^^]/2, respectively. Using Eqs. ([2|-(|6| 
we obtain, (see Supplementary Material for more details). 
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where an{0) = dEn{KO)/dkx\f.^Q and bn{0) = 
dEni)^^0) / dky\^^Q. Thus, the dependence of an{0) and 
Pnifi) on the magnetization orientation can be obtained 
by computing the k-space gradient (velocity) of the ab- 
initio energy bands in the vicinity of the F point. The 
functional forms of an{0) and bniO) conform to the sym- 
metry requirements (see Supplementary Material). 

Figure [3] shows the magnetization dependence of the 
Bychkov-Rashba-type and Dresselhaus-type SOC param- 
eters for the interface bands. The SOC parameters ex- 
hibit an oscillatory behavior as a function of the magne- 
tization orientation. The angular dependence of the SOC 
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FIG. 2: Spin-orbit-field "butterflies". Calculated spin-orbit 
fields for the magnetization along [110] (left) and [110] (right). 
The polar plots of the spin-orbit coupling strength (w/k) in 
the units of eVA as well as the corresponding vector fields 
w(k) are shown for the band n = 1 and the momentum con- 
tours of (a) k = 7r/100(i; (b) k = n/Sd; (c) k = 7r/5d. The 
lengths of the direction vectors have been rescaled. 



parameters is stronger for band n = 1 than for n = 2. 
In particular, for the case of band n = 1 the Bychkov- 
Rashba-type SOC parameter can even change its sign 
when the magnetization is rotated in the plane. This 
leads to the sign change of the product ai^i when the 
magnetization is rotated from [110] to [110] and produces 
the flipping of the SOF symmetry axes shown in Fig.|2|L). 
For band n = 2 the angular dependence is weaker, the 
product 0^2/^2 does not change its sign and the symmetry 
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axis of the SOF is preserved, being independent of the 
magnetization orientation. 

When considering the dependence on the transverse 
electric field, the behavior is opposite. Indeed, while the 
SOC parameters corresponding to band n = 1 change 
very little with for band n = 2 the changes in the 
magnitudes of 0^2 and f32 are appreciable. This disparate 
behavior is a consequence of the different nature of these 
two bands. Band n = 1 originates mostly from the two 
Fe atoms in the interface unit cell and, therefore, its cor- 
responding SOF is more sensible to the changes in the 
magnetization direction. However, the electrostatic con- 
trol of the SOC parameters is dominated by the elec- 
tric field influence on the pd bonding between As and Fe 
atoms. Consequently, the SOF corresponding to band 
n = 2, which comes mostly from the interfacial As atom, 
its neighboring Ga and the Fe atom above Ga, exhibits 
a stronger dependence on the electric field. 
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FIG. 3: Calculated magnetization and electric field depen- 
dence of the spin-orbit coupling parameters. The Bychkov- 
Rashba-type an and the Dresselhaus-type /3n spin-orbit pa- 
rameters for the interface bands n = 1, 2 are shown as a func- 
tion of the in-plane magnetization orientation and for different 
electric fields. 



sponding to the As-surface bands (n = 3, 4), which due 
to their surface nature possess stronger SOFs. 
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4(+) 


^(+) 






1 


-0.42 


-6.26 


-11.32 


4.32 


2 


-42.51 


1.82 


-57.94 


-1.51 


3 


-620.24 


-88.74 


-597.56 


-89.62 


4 


680.09 


95.61 


697.58 


103.15 



TABLE I: Band-resolved expansion coefficients of the 
Bychkov-Rashba-type and Dresselhaus-type spin-orbit cou- 
pling parameters in meV A units. The parameters are in the 
range of what is found in semiconductors [2]), 

If the magnetization is perpendicular to the interface 
plane, the first order correction to the energy vanishes 
and the methodology used for the case of in-plane mag- 
netization does not apply. To extract useful information 
about SOFs we first note that the in-plane components 
of the spin appear due to SOFs only (without spin-orbit 
coupling the electron spins would be fully polarized in 
the growth direction). Since the exchange field domi- 
nates over SOC, the spin is still quantized largely along 
the magnetization direction, so the expectation values of 
the transverse components of the spin {s)n correspond- 
ing to the nth band can be obtained by considering Hso 
as a perturbation. First order perturbation theory gives 
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where Axe is the exchange splitting energy and h = 1. 
Using these approcimate relations, we can determine the 
pattern of but not its magnitude. Figure[4]shows w{k.) 
(right parts), and its rescaled magnitude w = |i^(k)| (left 
parts) in the units of the exchange splitting Axe foi" the 
interface band n = 1. The fields have been computed on 
two contours, k = 7r/25d (a) and k = ir/Sd (b). Similar 
to the in-plane magnetization case, when the magneti- 
zation is perpendicular to the layers, the SOF close to 
the r point resembles the interference of the Bychkov- 
Rashba-type and Dresselhaus-type SOCs. The effect of a 
transverse electric field quantitatively modifies the SOF 
and is more pronounced for larger k values [see Fig. [4]: 



In Table I we list the expansion coefficients of the SOC 
parameters [see Eqs. (23-24) in the Supplementary Ma- 
terial] , 



an - +5i+)cos(2^), 
pn ^ 4"^+5i-)cos(2^), 



(10) 
(11) 



for zero electric field. From A\^^ ^ one extracts the 
magnetization-independent part, whereas the B^^~'^ pa- 
rameters control the leading contribution (higher order 
coefficients are about two orders smaller) to the angular 
(magnetization orientation) dependence of the spin-orbit 
parameters. In addition to the interface bands (n = 1, 2) 
we have also included the expansion coefficients corre- 



To summarize, we introduced a method to calculate 
spin-orbit fields from first principles and applied it to 
the Fe/GaAs interface. We found the the spin-orbit fields 
depend strongly on the Fe magnetization direction. This 
finding should be particularly important for lateral and 
tunneling magnetotransport anisotropics of ferromagnet- 
nonmagnet slabs. 

We thank D. Weiss, G. Bayreuther, C. Back, 
G. Woltersdorf for useful hints related to experimental 
ramifications of the presented theoretical concepts and 
F. Freimuth, Y. Mokrousov, G. Bihlmayer, J. Spitaler 
and P. Novak for helpful discussions regarding the cal- 
culations. This work has been supported by DFG SFB 
689. 



5 



(a) 



[110] 




E = -l V/nm I 
E = {) 

E = 1 V/nm 







(b) 



0.06 
0.03 


0.03 
0.06 



[110] 








FIG. 4: Electric field control of the spin-orbit coupling for 
magnetization perpendicular to the plane. Polar plots of the 
spin-orbit field strength (left panel) corresponding to the in- 
terface band n = 1 on the contours (a) k — 7r/25(i and 
(b) k — Tr/Sd, for transverse electric fields -1, 0, 1 V/nm. 
The momentum-dependent expectation values of the in-plane 
spin vectors {s) on the two corresponding /c-contours is shown 
in the right panel for zero electric field. The size of the vectors 
has been rescaled. 
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Magnetic control of spin-orbit fields: a first-principles study of Fe/GaAs junctions 
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I. GENERAL THEORY. 

We describe our symmetry-based method to obtain the in-plane spin-orbit fields (SOF) from first-principles data. 

In the presence of spin-orbit coupling (SOC) the energy bands become anisotropic with respect to the magnetization 
orientation m = (cos^, sin^, 0). The dispersion relation of the nth energy band can, in general, be written as 

£;,(k, 0) = TniK 0) + ^n{0) + A£;^°(k, o), (i) 

where T^(k, ^) is even in k. The k-independent term A^(^) describes the overall energy shift at the F point and is 
related to both the exchange energy splitting and the intrinsic atomic SOC. The SOC induced by lack of inversion 
symmetry generates the contribution A£'^°(k, ^), which is odd in k. By time reversal symmetry, the energy bands 
must remain invariant under the simultaneous inversion of the wave vector and the magnetization direction, i. e., 
En{k, 0) = k, 6 -\-7r). This, together with Eq. (1), allows for expressing the SOC correction in terms of the band 
energy, 

AE-(k,^)=^"('-'^^ (2) 

In our system not only the interface confinement but also the magnetism reduce the symmetry and lift the band 
degeneracies at the F point, as shown in Fig. 1 of the article.^ This allows for a description of the SOC in terms of 
Pauli matrices. As discussed in the article, the most general SOC Hamiltonian that exhibits C2v symmetry can be 
written as, 

nso = Wn{KO)-cr, (3) 

where cr are Pauli matrices and Wn (k, 6) , whose components are 

'^nxO^i^) — Vn{kxi ky, 0)ky ; Wny (k, 6) = fin{kx,ky, 0) kx , (4) 

is the spin-orbit field (SOF) vector. The x and y axes refer to the principal symmetry axes of the C2v symmetry. In 
our case x = [110] and y = [110]. 

The most general functional form of r]n{kx^ ky, 6) and fin{kx, ky, 6) which is compatible with the C2v symmetry of 
the SOF is given in Eq. (2) of the article. The functions r]n(kx,ky,6) and firi{kx^ky^6) are even in both kx and ky. 
Therefore, the SOF components must obey the symmetry relations, 

'^nx{kxi ky^ 6) = Wrtx{~kxT ky^ 0) ; Wny{kx-)ky^O) = —Wny{—kxiky^O). (5) 

Since the exchange field dominates over SOC, the contribution of Hso to the energy can be treated within first 
order perturbation theory, which yields the energy correction 

A£;;° (k, 6) = (J [wnx (k, 6) cos 6 + Wny (k, 6) sin 0] , (6) 

where a indicates whether the spin of the unperturbed band (i. e., in the absence of SOC) is parallel (cr = 1) or 
antiparallel (cr = — 1) to the magnetization. Away from band anticrossings the effect of SOC on the spin orientation 
is small and does not change the spin character of the band. Therefore, away from anticrossings, we still can reliably 
identify a = 1 (cr = — 1) with the spin- up (spin-down) character of the perturbed energy band (i. e., in the presence 
of SOC). 

Equation (2) together with Eq. (6) connect the SOF with the band energies. However, in order to extract Wnx{K 0) 
and Wnyi}^^0) separately, additional symmetry constraints are needed. The Fe/GaAs interface exhibits a structural 
C2v symmetry with the refiection planes (110) and (110) and [001] as the two-fold rotation axis (see Fig. 1). But 
when the magnetism of Fe is included, the symmetry is reduced and the C2v symmetry of the energy contours 
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FIG. 1. Schematic of the atomic structure of an As-terminated Fe/GaAs interface. The symmetry of the interface is that of 
the point group C2v, containing the twofold rotation axis C2 parahel to the growth direction [001] and the two mirror planes 
(lIO) and (110). 

in the momentum space is destroyed. Nevertheless, additional symmetry constraints on Wnxi^^^) ^"^^ '^ny(k, 6') 
are found by considering the magnetization as a symmetry variable. Note, for example, that (we consider only 
< ^ < 7r/2) the spatial reflection with respect to the plane (110) is equivalent to the rotation of the magnetization 
by —20 around the [001] axis (see Fig. 2).^ This implies that the energy contours in k-space must fulfill the relation 
En{—kx,kyO) = En{kx,ky, —6), which according to Eqs. (1), (2), and (6) yields. 



'^nx{kx-)ky^ 6^ — Wfix{ kx-iky^O^ 5 '^ny{kx-)ky^ 0^ — '^nyi kx^ky^O^. 

By combining Eqs. (5) and (7) we obtain the following symmetry constraints on the SOF components, 

"i^nxikxi ky^ 6) = Wnx{kxi ky, —6) ; Wny{kx-)ky^6) = Wny{kx-)ky^ —6\ 
which together with Eq. (6) leads to 



Wnxi}^, 0) = a 



A£;-(k,^) + A£;-(k,-^) 



2cos^ 



; Wny{KO) = (J 



A£;-(k,^)-A£;-(k,- 



2 sin 6* 



With the help of Eq. (2) we can rewrite the above relations as, 

'En{K 0) - ^n(-k, 0) + En{-kx, ky, 0) - En{kx, -ky, 0) 



and 



'^nxikx-) ky-i 0^ — (7 



Wny{kx,ky,6) = a 



4cos^ 



En{K 0) - £;n(-k, 0) - En{-kx, ky, 0) + En{kx, -ky, 0) 



4 sin 6* 



(7) 
(8) 

(9) 
(10) 
(11) 



which allows us to extract the SOF directly from ab-initio data.^ 

In the particular cases 6 = 7r/2 and ^ = the numerators and denominators in Eqs. (10) and (11) vanish and the 
SOF can be obtained from L'Hopital's rule. 



II. LINEAR SPIN-ORBIT FIELDS. 



In the vicinity of the F point the SOF is linear in the wave vector components kx and ky. According to Eqs. (2) 
and (3) of the article, the linear contributions to the SOF are given by 



(k, 0) = ,7(0) {e)ky ■ (k, 9) = {e)K, 



where 



dw„,{k,e) 



dw^yi'k, e) 



(12) 



(13) 
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FIG. 2. (a) Top view of the investigated Fe/GaAs interface as seen from the interface plane of As atoms. The red spheres 
represent four Fe atoms nearest to the As atom (yellow sphere). The arrows represent the Fe magnetic moments aligned along 
the magnetization direction m = (cos 0, sin 0, 0). (b) The structure of (a) after a reflection with respect to the mirror plane (110) 
(recall that spins and magnetic moments are axial vectors). Alternatively, one can also obtain (b) from (a) by a rotation of the 
magnetization by the angle —26 around the [001] axis. This leads to the symmetry relation En( — kx, kyO) = En(kx, ky, —9). 



Taking into account Eqs. (10) and (11) we obtain, 



bn{0) 



where 



dkr 



bn{0) : 



sinO 



(14) 



(15) 



Thus, the size of the linear SOC induced by inversion asymmetry can be estimated by computing the k-space gradient 
(velocity) of the ab initio energy bands in the vicinity of the F point. 

The general form of coefficients an{0) and bn{0) can be deduced from symmetry considerations. By time reversal 
symmetry, £'n(k, 6) = En{—\^, -\- tt) and, therefore, the following relations 



an{e) = -an{e + tt) ; bn{e) = -bn{e + 

must hold. On the other hand the constraint in Eq. (8) leads to 

an{0) = -an{-e) ■ bn{0)=bn{-e). 



(16) 



(17) 



Furthermore, since both an{0) and bn{0) are periodic functions of 0, they can be expanded in Fourier series. Thus, 
under the symmetry constrains in Eqs. (16) and (17) the corresponding Fourier expansions for an{0) and bn{0) reduce 
to 



{d) = J2(inl sm[(2/ + l)^] 



(18) 



and 



b„{e) = Y,bnicos[{2l + l)0], 



(19) 



respectively. For the bands investigated in the article we have found that for accurately fitting the band energy 
gradient at the F point it is enough to consider the first three terms in the series expansions of an and bn- Thus, one 
can further approximate Eqs. (18) and (19) as 



and 



o^n (^) = CLnO sm 6 + tt^i sin(3^) + a^2 sin(5^) 
bn{^) = bno COS + bni cos(3^) + 6^2 cos(5^). 



(20) 
(21) 
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FIG. 3. Calculated spin-orbit coupling coefficients (a^ and bn) for the Fe/GaAs slab corresponding to the interface band n = 2 
as a function of the magnetization orientation. The lines are fittings of Eqs. (20) and (21) to the ab-initio data (dots). The 
data are obtained by using Eq. (15). Green dots and dotted (green) lines represent the case of zero electric field, while red 
circles and sohd (red) hues correspond to the electric field E = —1 V/nm, and blue squares and dashed (blue) hues to the field 
E = 1 V/nm . The magnetization orientation is measured with respect to the [110] crystallographic direction of GaAs. 



respectively. 

For illustration we show in Fig. 3 a comparison between the magnetization direction dependence of the n = 2 
band-energy gradient at the F point obtained from the ab-initio data and the corresponding fittings using Eqs. (20) 
and (21). The good agreement corroborates that the higher order terms in the series expansion can here be indeed 
neglected. The fittings are excellent even in the presence of a finite electrostatic field along the [001] axis, as can 
be appreciated from Fig. 3. This is because the presence of such an electrostatic field does not affect the symmetry 
considerations used in the derivation of Eqs. (18) and (19). 

As discussed in the article the Bychkov-Rashba-type and Dresselhaus-type SOC parameters, an and /3n respectively, 
can be introduced through the relations r]n^ = f^n — and /j,n'' = l^n + Cin- Therefore, according to Eq. (14), we 
obtain. 



an (0) COS — bn (0) sin 



sin(26') 



an [0) cos -\-bn {0) sin 6 



sin(26') 



(22) 



which allow us to extract the dependence of an{0) and /3n{0) on the magnetization orientation from the ab-initio 
data. By combining Eqs. (20), (21), and (22) we obtain the following expressions for the angular dependence of the 
Bychkov-Rashba-type and Dresselhaus-type SOC parameters. 



an{e) = a + 5^+) cos(2^) + cos(4^)] , 

t5n{0) = a [4-) + cos(2^) + Ci-) cos(4^)] , 



with and ci^^ given by 



(23) 
(24) 

(25) 



B^^^ = a„i + a„2 ± (b„2 - 6ni), 



(26) 
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and 

Ci,^^ =an2Tbn2, (27) 

respectively. 

In the bottom panel of Fig. 3 of the article we show the magnetization direction dependence of the SOC parameters 
for the case of the interface band n = 2 and three different values of the electrostatic field across the heterostructure. 
For the plots we used Eqs. (23)- (27) with the expansion coefficients ani and bni (/ = 0, 1, 2) obtained after the fitting 
in Fig. (3) and took into account that band n = 2 exhibits spin-down character near the F point (i. e. a = —1). 
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^ Similarly, for7r/2<^<7ra reflection with respect to the mirror plane (110) is equivalent to a rotation of the magnetization 
of the initial system by an angle n — 20 around the [001] axis. In such a case one obtains En(kx, —kyO) = En(kx, ky,7r — 6), 
which is the time reversal of the situation discussed in the text. 

^ A similar analysis for the case 7r/2 < < tt leads to the same expressions [Eqs. (10) and (11)] for the SOFs. Therefore, 
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